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Abstract 

We study the problem of option replication under constant proportional transaction costs 
in models where stochastic volatility and jumps are combined to capture market’s important 
features. In particular, transaction costs can be approximately compensated applying the 
Leland adjusting volatility principle and asymptotic property of the hedging error due to 
discrete readjustments is characterized. We show that jump risk is approximately eliminated 
and the results established in continuous diffusion models are recovered. The study also 
confirms that for constant trading cost rate, the results established by [14] and by [25] are valid 
in jump-diffusion models with deterministic volatility using the classical Leland parameter. 
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1 Introduction 

Many suggested mathematical models for stock prices have been trying to capture important mar¬ 
kets features, e.g. leptokurtic feature, volatility clustering effect, implied volatility smile. These 
market properties are tractable in stochastic volatility models. However, it is worth noticing 
that diffusion-based stochastic volatility models, where the market volatility can fluctuate au¬ 
tonomously, can not change suddenly and as a result, they could not take into account sudden and 
unpredictable market changes. Hence, for a realistic setting, the continuity assumption of stock 
price should be relaxed. In fact, as discussed in [27], the presence of jumps in the asset price 
can be recognized as the presence of participants in the option market. In order to construct an 
realistic extension of the famous Black-Scholes framework, it is reasonable to suppose that good 
or bad news arrive according to a Poisson process. The changes of asset price are described by the 
jump-sizes and between two jump times, the asset price follows a geometric Brownian motion as in 
the classical Black-Scholes models. Such a combination is called a jump-diffusion model. As shown 
in [1G], jump-diffusion models not only fit the data better than the classical geometric Brownian 
motion, but also well reproduce the leptokurtic feature of return distributions. See [27, 26, 16] and 
the references therein for detailed discussions. 
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It is well-known that in complete diffusion models, options can be completely replicated using 
the delta strategy adjusted continuously. However, it is not the case for models with jumps. In 
fact, jumps risk can not be released completely even using continuous time strategies and the only 
way to hedge perfectly an option against jumps is to buy and hold the underlying asset. In other 
words, the conception of replication does not indicate a right framework for risk management and 
hedging in the presence of jumps as in diffusion-based complete market models where Black-Scholes 
theory plays a central role. Therefore, discrete hedging for jump models is practically important 
and asymptotic properties of the hedging error are not a trivial task to show. 

The situation becomes more challenging if one takes into account transaction costs which are 
needed for hedging activities in practice. Such a consideration is realistic and has been attractive to 
researchers for last years. In the absence of jumps, [17] proposed a simple method to compensate 
trading costs, which is a modification of the well-known Black-Scholes PDE where volatility is 
artificially increased as 

a 2 = a 2 + ann 1 / 2 - 01 ^/8/n, (1) 

where n is the number of revisions and Kn~ a , 0<ct<l/2is the cost rate. He claimed that for 
a = 0 (constant rate) or a = 1/2, using the corresponding discrete delta strategy can produce 
a portfolio whose terminal values V" approach to the option payoff h(S i) as revision becomes 
frequent. 

However, [14] showed later that the Leland statement for constant transaction cost is not 
mathematically correct and hedging error in fact converges to a non-zero limit min(Si, K ) — J(S i) 
as the portfolio is frequently revised, where J(S'i) is the limit of cumulative costs. The rate of 
convergence was then investigated by [25] including a characterization of asymptotic distribution 
of the corrected replication error. In particular, it was shown that the sequence 

n 1 / 4 (H" - h(Si) - min (S U K) + kJ(Si)) 

weakly converges to a mixed Gaussian variable as n —» oo. This result has initiated many further 
studies for the case of constant cost in different directions: studying the problem for option with 
general payoffs, using non-uniform readjustments to accelerate the rate of convergence by [18, 19, 5], 
considering the problem in more general models e.g. local volatility by [?], trading costs based on 
the traded number of asset by [?]. Recently, [23] have studied the problem in stochastic volatility 
frameworks using a simpler form for adjusted volatility. It turns out that increasing volatility 
principle is still helpful for controlling losses caused by trading costs which are proportional to the 
trading volume in different general situations. Furthermore, a connection to asset hedging in high 
frequency markets, where the form of ask-bid price may be an essential factor for deciding laws of 
trading costs, is also discussed. We refer the reader to the above papers and the references therein 
for a full discussion on Leland’s approach. 

To the best of our knowledge, there is few study about the trading costs for models with jumps 
following the Leland spirit in discrete time setting. In fact, it is very intuitive to think that in the 
presence of transaction costs and/or jumps in the asset price, the option is more risky and should 
be evaluated at higher price than that in the absence of these risks. However, a more expensive 
option price would be equivalent to say that there has been some increase in its volatility values. 
That is the essential intuition behind the Leland algorithm. 

The aim of the present note is trying to build an a bridge from the existing results in continuous 
diffusion models studied by [ L4, 25, 18, 19, 5, 23, 24]... to discontinuous models where jumps are 
allowed in asset price and/or volatility. In fact, we try to capture not only the dependence structure 
(using stochastic volatility) but also short term behaviors of the stock price due to sudden market 
changes as well 1 . This combination is expected to model a market general enough for some practice 
purposes since stochastic volatility models well complement models with jumps [16]. 


1 This fact partially explains why jump-diffusion models are, in general, considered as a good choice, especially 

in short-term situations. 
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Figure 1: The corretor of hedging error on the left and limit of trading costs on the right with 
strike K = 5 


It is well-known that jump risk in hedging problem is challenging to handle can not be released 
completely even in simple framework e.g jump-diffusion models and continuous adjustments are 
possible. In this paper, we show that impact of jumps can be partially negligible under some mild 
condition on jump sizes. In fact, asymptotic distribution of the hedging error is independent of 
jumps and consistent with those established in [23]. The same thing is true when jumps are allowed 
in both asset price and its volatility. Such general frameworks provide the ability to explain large 
movements in volatility, which happen during crisis periods [7, 8]. Moreover, it is important to 
emphasize that the present study also proves that for constant cost rate, the Kabanov-Safarian- 
Pergamenshchikov results in [14, 25] are true for jump-diffusion settings. This is also a new 
contribution to the literature. In fact, it is quite straightforward to extend this results for jump- 
diffusion models to the case when option payoff is defined by a general convex function holding 
some decaying assumption studied in [ 3, 19, 5] and in other works. 

The remainder of this paper is organized as follows. We shortly present the key idea behind 
the Leland suggestion in Section 2. Section 3 is devoted to formulate the model and present our 
main results. General stochastic volatility models with jumps will be discussed in Section 4. We 
discuss the special case when volatility is constant in Section 5 and present a numerical example in 
Section 6. Proof of main results are reported in Section 7 and some useful Lemmas can be found 
in the Appendix. 


2 Which form for adjusted volatility in Leland’s algorithm? 


To explain the key idea in the Leland’s algorithm we assume that the stock price is given by 
d St = aStdWt and the interest rate is zero so that S' is a martingale under the risk-neutral measure. 
Under the presence of proportional transaction costs, it was proposed by [17] then generalized by 
[14, 15] that the volatility should be adjusted as in equation (1) in order to create an artificial 
increase in the option price to compensate possible trading fees. This form is inspired from the 
observation that the trading cost S t \C x (t t , S t _ ) — C x (t t , S t )| interval of time [t i _ 1 ,ti] can 
be approximated by 


K nS t ._ 1 C xx (t ti _ i , S t ._ i )\AS t . \ ss Kn aSl_C xx (t u _ i ,S u _ i )E\AW u 


( 2 ) 
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For simplicity, we assume that the portfolio is revised at uniform grid ti = i/n , i = 1,..., n of the 
option life interval [0,1]. Taking into account that E|AW / t ./(At i ) 1 / 2 | = -^/2/7r one approximates 
the last term in (2) by Kna^JZ/TtlyAtiY^S 2 which is the cost paid for portfolio 

readjustment in [t i _ 1 ,ti\. Hence, by the standard argument of Black-Scholes theory, the option 
price inclusive of trading cost should verify 


Ct(t ti _ v S u 


JA ti + -<7 2 £j 
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U- 


C xx(t ti _ i ,S u _ i )At i + K n CT V /277r(Atj) 1 / 2 5' 2 _ i C' ra (t t ._ i ,5' t ._ i ) = 0. 


Since At, = 1/n, one deduces that 


+ \ (CT 2 + K„cr v / n8/V) St i _ 1 C xx (t u _ 1 ,S t ._ 1 ) =0, 


which implies that the option price inclusive trading cost should be evaluated by the version of the 
BS PDE 

C t (t, x ) + \d 2 x 2 C xx (t , x) = 0, (7(1, x) = h(x), (3) 


where the adjusted volatility a defined by (1). This Leland’s observation is important for practical 
purposes since it requires only a small change in parameter of the well-known Black-Scholes frame¬ 
work that is wildly used in practice. However, when volatility is random, say a = cr(yt) where y t is 
another random process, the strategy is no longer available from the Cauchy problem (3). In fact, 
our model is now called stochastic volatility (SV) models. Pricing and hedging in such contexts are 
intrinsically different from the Leland model and even in the absence of trading costs, the option 
price strongly depends on volatility level and future information of the volatility process y. To see 
this issue, let us assume that y is a Markov process. Then by iterated expectation the option price 
exclusive of trading costs is given by 


C(t,x,y) = E,[C BS (t,x- : K,a\y t = y)\, (4) 

where K is the strike price and er is the averaged volatility defined by 

1 f 1 

v 2 = J (^) ds - 

It means that option price is the average of Black-Scholes price on all possible future trajectories of 
the volatility process y, which in reality can be not observed directly. Thus option price and hence 
hedging strategy is very complicated to get in a closed form and usually studied via asymptotic 
analysis. 

Hence, in the presence of transaction costs, the above discussion emphasizes that the well- 
known form (1) for adjusted volatility in Leland’s algorithm is no longer helpful from a practical 
point of view. Furthermore, even for local volatility models when a depends on time and the spot 
price 5, it is technically difficult to show the existence of a solution to (3) since the operator is not 
uniformly parabolic in such cases. In addition, estimates for derivatives of option price, which are 
essential for approximation analysis, are not an easy task, see [18]. 

Fortunately, a deep study on the approximation of the hedging error shows that the limit of the 
replication error does not strongly depend on the form of adjusted volatility but only on the last 
term K n ^/n whenever the latter product diverges to infinity. This important observation means 
that a simpler form gn n ^/n with some constant g, can used to obtain the same asymptotic property 
for the hedging error. This modification is fully investigated for SV models in [23]. Note that the 
option price and hedging strategy are easily obtained if the new form is applied for the Cauchy 
problem (3) while the classical one (1) is impossible in practice. In this paper we show that this 
suggestion still useful even when jumps in asset price and/or in stochastic volatility are taken into 
account. 
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We conclude the section by mentioning some important features of hedging with jumps models. 
First, jump risks are not covered by simply using the classical delta strategy even with a continuous 
time policy. In other words, in the presence of jumps delta hedging is no longer optimal. Second, if 
jumps are allowed in the stock price then one should distinguish two types of hedging errors: one 
is due to market incompleteness concerning jumps and the other one is due to the discrete nature 
of the hedging portfolio. These two types of hedging errors have different behaviors. 

The literature of discrete hedging with jumps is vast and we only mention to [29, 30] for recent 
achievements. Seemingly, none of these mentioned papers discussed about trading costs. 


3 Model and main results 


3.1 The market model 

Let ^17, F \, (7i)o< t <i, be the standard filtered probability space with two standard independent 

(J r t) 0<t <i adapted Wiener processes (' Wj 1 ^) and (Ik/ 2 ' 1 ) taking values in R. Consider a financial 
market consisting of one non-risky asset set as the numeraire and the risky one (e.g. stock ) S t 
being assumed to take jumps at random times (tj) , > . | . Relative changes in value of the risky 
process at the jump time r 7 are characterized by the sequence of i.i.d. random variables (£j)j>i- 
The jumps-size process is then defined by 


— Col{o} + (t) > Co — 0. 


J>i 


(5) 


We assume further that between the jumps times, St follows the classical geometric Brownian 
motion but with stochastic volatility. More precisely, we suppose that the risky asset dynamics is 
driven by the following equations 


dS t = S t - (b t dt + a ( y t ) dIF t (1) + dCt) , 
d y t = ai ( t , y t ) d t + a 2 (t, y t ) d Wp\ 


( 6 ) 


where (t. = C? is the jump part characterized by the Poisson process N t with intensity 

parameter 9 and S t - = lim stt S s . We assume that the coefficients at,i = 1,2 are locally Lipshitz 
and linearly growth functions, which provide the existence of the unique strong solution y to the 
second equation [10]. 

In this paper, all sources of jump randomness mentioned above: two Brownian motions, the 
Poisson process N t and the jumps sizes (Cj)j>l are assumed to be independent. Note that limiting 
results of replication error will change only on its asymptotic distribution if the two Brownian 
motions are correlated, see [2 ]. 

Let us explain how the asset fluctuates. In fact, it can be observed that the risky asset price 
St changes continuously most of the time, but large jumps may occur from time to time. This 
important fact can not be adequately taken into account in the classical Black-Scholes context 
or other pure diffusion-type models. Moreover, it is easy to see that system (6) has the unique 
solution 


i rt rt i ft N t 

/ b t dt+ er(y s )dIF s (1) — - / a 2 (y s ) ds + ^ In (1 + £,•) 
Jo Jo Z Jo j =1 


(7) 


which means that the stock price is simply a product of the geometric Brownian motion with 
stochastic volatility Sq exp | J 0 * a ( y s ) dW^ — \ /J cr 2 (y s ) ds + / 0 * 6*ds| and the independent jump 


N t 


part I] (1 + £j)■ The jump sizes should satisfy the natural condition 
i=i 

> —1, for all i = 1, 2,_ 


( 8 ) 
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to guarantee the positivity of the stock price. We easily observe that at the jump time r 7 the 
relative change in value of S is given by A S Tj = S T - (£,■ + 1). Let J (d t x dz) be the random 

Poisson measure generated by the compound Poisson process Q, that is 

J (M X A) = ^ 1 { Ti <t, £j€A}- (9) 

j> 1 

It is well known that J (df x dz) is a er-finite jump measure whose Levy measure (intensity measure) 
is defined by v = OF (dz) df, where F(-) is the common distribution of jump sizes (£,). The Levy 
measure can be interpreted as the average number of jumps per unit of time. For convenience, we 
use the notation v(dz) = 0F(dz) for its density and denote the compensated stochastic Poisson 
measure by J (dt x clz) = J (dt x dz) — v(dz)dt. The jump measure J (di x dz) permits to define 
the stochastic Poisson integral (integral with jumps) 

n N t 

f(s,z)J(dsxdz) = ^2f(Tj,Zj), ( 10 ) 

i - =1 

for any predictable process / Adapting from the pricing principle we assume that the drift 

b = — 0E£i so that the stock price is now a local martingale given by 

d S t = S t - (a ( y t ) dW t (1) + J zdJ(dt x dz)) . (11) 

Remark 1. The drift b t plays no role in approximation. In fact, our asymptotic results are valid 
for any b t holding sup 0<t<1 E|b t | < oo. 


Volatitlity paths 



Asset price paths 



Figure 2: Asset price and its volatility paths in Hull-White model with Gaussian jumps in asset. 


Remark 2. In this paper we do not discuss about the problem of change of measure and jump risk 
but accept the free-risk assumption of asset dynamics as the starting point. Clearly, jump-diffusion 
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suggestion leads to an incomplete market, which is also an important feature of stochastic volatility 
settings. Hence, there are many ways to choose the pricing measure throughout the Girsanov 
technique. Such a procedure makes an essential change not only on the diffusion but also on the 
jump part of the asset dynamics [27, 22]. In [16], the rational expectations equilibrium is used 
to obtain a simple transform from the original physical probability to a risk-neutral probability so 
that many assets (bonds, stocks, derivatives on stocks) can be simultaneously priced in the same 
framework. 


3.2 Assumptions and examples 

The following condition on jump sizes is accepted in our consideration: 
(Ci) The common distribution F of jump sizes satisfies 


[ z 2 F(dz) < oo and [ —-— F(dz) 
Jr J-i 1 + 2 


< oo. 


The first integrability condition is nothing than the condition of finite variance for jump size dis¬ 
tribution while the second one is equivalent to E(l + £) -1 < oo. These are quite weak constraints 
and automatically fulfilled in the Merton jump-diffusion model [22] where jump size distribution is 
assumed to be log-normal. In [16], within an equilibrium-based setting, log-exponential distribu¬ 
tions are suggested for jump sizes to obtain the convenient feature in analytical calculation. Again, 
this family of jump size distributions verifies condition (Ci). 

Let us turn our attention to volatility assumption. Following [23] , we assume that the volatility 
process satisfies the integrable condition. 

(C 2 ) The function a is twice continuously differentiable such that 

0 < cr min < inf <j(y) and sup E m&x{a(y t ), \<j\yt)\} < 00 . 
y£ R o<t<i 


In fact, condition (C 2 ) is fulfilled for almost of the famous stochastic volatility models, see [23] 
for more discussions. 


Remark 3. It is important to note that in the present setting, the combination of stochastic 
volatility and jumps means that the asset price is not a Levy process but a semi-martingale. As 
mentioned in [23], finite moments of the asset process in a general dynamics of volatility process 
are not insured, see [1, 20]. This crucial feature prevent us from making a L 2 -based approximation 
as in deterministic volatility models [If, 18, 19]. Therefore, for approximation procedures, we 
follow the approach in [23, 24] based on a truncation technique and convergence results obtained 
are guaranteed in probability. 

We conclude this subsection with some well-known stochastic volatility models with jumps, see 
[27] and Section 4 for more examples. 

The Bates models: The Bates models is a jump-diffusion stochastic volatility models obtained 
by adding proportional log-normal jumps to the Heston stochastic volatility model: 

dS t = Sffpdt + y/yidWf + dZ t ); dy t = a(m - y t )dt + bffytdWf, (12) 

where W s , W v are Brownian motions with correlation p and Z is a compound Poisson process with 
intensity 0 and log-normal distribution. Condition (Ci) is clearly verified since jumps follow the 
log-normal law. Bates’s models exhibit some very nice properties from a practical point of view. 
Firstly, the characteristic function of the log-price is available in a closed-form, which is important 
for pricing purposes. Secondly, the implied volatility pattern for long term and short term options 
can be adjusted separately [27]. 
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Ornstein-Uhlenbeck stochastic volatility models: It is possible to introduce a jump compo¬ 
nent in both price and volatility processes. Such models are suggested by Barndorff-Nielsen and 
Shephard to take into account leverage effect: 

S t = S 0 exp X t \ dX t = (p + /3af )dt + cr t dW t + pdZ t \ dcr \ = —0of dt + d Z t . (13) 

If p = 0 the volatility moves with jumps but the price process has continuous paths. The case 
p ^ 0, representing a strong correlation between volatility and price, provides the model flexibility 
but computation is now challenging. Remark that in this case a is not the ’’true” volatility since 
the returns are also affected by changes of the Levy process Z t . If jumps still follow log-normal 
law then condition (Ci) is fulfilled. 

3.3 Approximate hedging with transaction costs : main result 

We study the problem of discrete hedging in friction contexts using the increasing volatility prin¬ 
ciple as in Leland’s algorithm. A brief review on this literature can be found in [23]. See more in 
[25, 15, 19] for related results. In the present framework, we follow the setting in [23]. More pre¬ 
cisely, we suppose that for each successful trade, traders are charged by a cost that is proportional 
to the trading volume with the cost coefficient k. Here n is a positive constant defined by market 
moderators. Let us suppose that the investor plans to revise his portfolio at dates (U) defined by 

ti=g(i/n), g(t) = 1 - (1 - tf , (14) 

where n is the number of revisions. The parameter p is used control the rate of convergence of 
the replication error. The bigger value of p the more frequent near by the expiry the portfolio 
revisions. Clearly one gets the well-known uniform readjustment for p = 1. 


inhomogeneous trading times 



Figure 3: Trading times with different values of p 


To compensate transaction costs caused by hedging activities, the option seller is suggested to 
follow the Leland strategy defined by the piecewise process 
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(ti-i,S t - ^j (t ), 


(15) 
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where C is the solution to the following adjusted-volatility Black-Scholes PDE 

C t {t,x) + ^x 2 C xx (t,x) = 0 , 0 < t < 1; C(l,x) = h(x) := (x - K) + , (16) 

Here the adjusted volatility a 2 is given by 

S 2 {t) = Qyjnf'{t) with / = g~ l ■ (17) 

Motivation of this simple form is discussed in Section 2, see more in [23]. Now, using strategy 7 " 
requires a cumulative trading volume measured in dollar value given by T n = f S t | 7" — 7" |. 

Thanks to Ito’s lemma one represents the payoff as 

MSi) = C(0, So) + j d x (t, S t -)dS t + £ (d t (t, St) + l -a\y t )S 2 t .d xx {t , S t -)) d t 

+ J2 (c(t,s t )-d(t,s t -)-AStd x (t,s t -j), 

0<t<l 


where A St = St — S f - is the jump size of stock price at time t. The last sum of jumps can be 
represented as f* f R B (t, S t -, z) J (d t x dz), where 


B(t, x, z) = C (■ t , x(z + 1)) — C (t, x) — zxC x (t, x) (18) 

with the jump measure J(dt x dt) defined by (10). Define then I :i n = J 1 f R B (t, S t -, z) J (d t x dz). 
Assuming that the initial capital (option price) is given by Vq = (7(0, So) and using (16), one 
represents 2 the hedging error as 

vr - Hsj = \i hn + / 2jn - / 3>n - «r„, (19) 

where h,n = fg fit - a2 iVt )) S ?Cxx(t’ S t) dt and J 2,n = fg ( 7 - C x (t, S t -)j d S t . 

The goal now is study asymptotic property of the replication error V™ — h(S ± ). To describe 
asymptotic properties, let us introduce the following functions 


v(A,x) = ln ^ A ^ + g(A,x) = ~ \ and ^(A) x) = ip (v(A, x)), (20) 


where is the standard normal density function. As shown below, the rate of convergence of our 
approximation will be controlled by the parameter /? defined by 


\- i3: =WZT) < 1’ for 1Sf,<2 - 


( 21 ) 


Before stating our main result, let us emphasize that using an enlarged volatility which diverges 
to infinity implies that asymptotic property of hedging error strongly depends on trading times 
near by the maturity. But remember that jumps are rare events and hence, jumps near by the 
expiry can be omitted with very small probability. Therefore, jumps in such contexts do not much 
affect asymptotic property of the hedging error as revisions become more shorter. This important 
fact proves that increasing volatility as in [23] is still helpful for models with jumps. The below 
theorems are in fact the achievement of [23] for continuous stochastic volatility models. 


2 Note that for Lebesgue and Ito’s integrals one can replace S f - by S t 
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Theorem 3.1. Under conditions (Ci) — (C 2 ), the sequence of 

n 0 (Vr - fc(Si) - min (S U K) + kT (S uyi , q)) 

converges to a centered-mixed Gaussian variable as n tends to infinity, where the positive function 
r is the limit of trading volume defined as 

r +00 

r (x,y,g) = x A _1/2 ^(A,x)E \a(y)g~ 1 Z+ q{X,x)\ dA, (22) 

do 

in which Z is a standard normal variable independent of Si,yi. 

The term q( A, x) in the limit of transaction costs can be removed using the modified Leland 
strategy, so-called Lepinette’s strategy: 

7 1 = E (c x (ti-i,S tT i ) - J*'- 1 C xt (t,S t -)dt^ l (t 4 _ llt 4 ](i). (23) 

Theorem 3.2. Suppose that Lepinette's strategy is used for option replication. Then, under 
conditions (Ci) — (C 2 ) the sequence of ( V™ — h(S\) — 77 mm (Si, K)~) converges to a centered- 
mixed Gaussian variable as n tends to infinity, where r) = 1 — na(yi) g^ 1 ^/8 /tt . 

3.4 Super-hedging and price reduction 

As proved in [23], a suitable choice of g can lead to super-replication. 

Proposition 3.1. Let conditions (Ci) — (C 2 ) hold and a be a twice continuously differentiable 
and bounded function. Then there exists g* > 0 such that lim 7l _>. (>0 Vf > h(Si) for any g > g*. 
This property is true for both of Leland’s strategy and Lepinette's strategy. 

From Black-Scholes’s formula one observes that both strategy 7 " and 7 " approach to the buy- 
and-hold one as n —> 00 . This means that option is now expensive from the buyer’s point of view. 
In fact it is close to the buy and hold price Sq. In [25, 23] a simple method is suggested to lower 
the option price following the quantile hedging spirit. Let us adapt the main idea in these works 
for the present setting. Since S x - = S x almost surely, we define 

S e = inf {a > 0 : T(a) > 1 — e} , (24) 

where T(a) = P ((1 — n) min(5 , 1 , K) > (1 — a)S 0 ). The quantity S e is called quantile price of the 
option at level £ and the difference (1 — S £ )Sq is the reduction amount of option price (initial cost 
for quantile hedging). Clearly, the smaller value of S E is, the cheaper the option is. We show that 
the option price is significantly reduced, compared with powers of parameter e. 

Proposition 3.2. LetS s be Leland price defined by (24) and assume that the jump sizes are almost 
surely non-negative, i.e. 

> 0, a.s. Mj e N (25) 

and °max = suPj/gR o-(y) < oo . Then, for any r > 0, 

lim(l — S E )e~ r = +00 . (26) 

£—>-0 

Proof. Observe that 0 < S e < 1 and S E tends to 1 as e —> 0. Set b = 1 — k. Then for sufficiently 
small e such that S s > a > 1 — bK/S 0 one has 

1 - e > P(min(5 1 , K) > (1 - a)S 0 ) = 1 - P(5 1 /5 0 < (1 - a)). 
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( 27 ) 


Therefore, e < P (S 1 /S 0 < (1 - o)(l - k) *) = P (P 1 (^)£ 1 (y) < z a ), 

£ < P (SJSo < (1 - «)(1 - z^)- 1 ) = P (Pi(^)fi(y) < z a ), 
where z a = (1 — a)(l — z«) _ 1 e AE ^ 1 and 

£ t (a) = exp a(y s )dW^ 1) ~ \ cr 2 (j/ s )ds| and P t (£) = JJ(1 + £,-). (28) 

By (25), P t (£) > 1 for all t £ [0,1], which implies that the probability in the right side of (27) is 
bounded by P (£\(y) < z a )- Therefore, £ < P (£i{y) < z a ) ■ At this point, the conclusion exactly 
follows from Proposition 4.2 in [ ] and the proof is completed. □ 

4 General stochastic volatility models with jumps 

4.1 Introduction 

Stochastic volatility with jumps in price (SVJP) models have been very popular in the option 
pricing literature since they provide flexibilities to capture important features of returns distri¬ 
bution. However, empirical studies [7, 8 ] show that they do not well reflex large movements in 
volatility assets during periods of market stress such as those in 1987, 1997, 2008. In other words, 
SVJP models are misspecihed for such purposes. The studies also suggest that it would be more 
reasonable to add an extra component into the volatility dynamics so that this new factor allows 
volatility to rapidly increase. Note that such expected effect can not be generated by only using 
jumps in returns (as in jump-diffusion models) nor by using diffusive stochastic volatility. In fact, 
jumps in returns can only create large movements but they do not have future impact on returns 
volatility. On the other hand, diffusive stochastic volatility driven by a Brownian motion only gen¬ 
erates small increase via sequences of small normal increments. Empirical analysis in important 
works [7, 8 ] shows that incorporating jumps in stochastic volatility can provide rapid changes in 
volatility. 

It is important to note that introducing jumps in volatility does not means an elimination of 
jumps in returns. Although jumps both in returns and volatility are rare, each of them plays 
an important part in generating crash-like movements. In crisis periods, jumps in returns and 
in volatility are more important factor than the diffusive stochastic volatility in producing large 
increases. We refer the reader to [7, 8 ] for more influential discussions about financial evidence for 
motivation of the use of jumps in volatility. 

In this section, we study the problem of option replication under transaction costs in a general 
SV models with jumps in return as well as in volatility, which is clearly a generalization of the 
setting in Section 3. In such contexts, jumps in volatility can be also ignored as those in asset 
price, i.e. the results obtained in Section 3 are recovered when jumps are allowed in both asset 
price and volatility. 

4.2 Specifications of SV models with jumps 

Assume that under the objective probability measure, the dynamics of stock prices S are assumed 
to be given by 

d S t = S t - (b t (y t )dt + a(y t )dW^ ] + dCf) , d y t = aq (t, y t )dt + a 2 (t, y t ) dW t (2) + d£ t y . (29) 

Here, (f = and C t = are two compound Poisson processes. For a general setting, 

two Poisson processes N[ and two sequence of jump sizes (£J), r £ {S,y} can be correlated. Let 
us give some possible specifications for jump components. 
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(i) Stochastic volatility model (SV): Clearly, this corresponds to the case when there is no jump 
in both asset price and volatility, i.e. = 0,V£. This basic SV model has been 

widely investigated in the literature. The problem of approximate hedging under proportional 
transaction costs is studied in [23, 24]. Roughly speaking, adding some extra component 
generated by a diffusion to the returns distribution of a classic Black-Scholes setting gives a 
SV model. 

(ii) Stochastic volatility with jumps in volatility (SVJV): By allowing jumps in volatility process 
y one can obtain an extension of SV models, i.e. Q s = 0,Vf but Q 0. In such cases, option 
pricing implications are in fact inherited from SV models. 

(iii) Stochastic volatility with jumps in price (SVJP): Assume now that £ t s y 0 but Q = 0. This 
case is studied in Section 3. 

(iv) Stochastic volatility with common jumps in price and volatility (SVCJ): Suppose that both 
asset price and its volatility in a SV model are influenced by the same extra random factor 
modeled by a compound Poisson process. In other words, jumps in asset price and in volatility 
are driven by the same compound Poisson process (jf = Cjf. 

(v) Stochastic volatility with state-dependent and correlated jumps (SVJJ): This is the most gen¬ 
eral case for the present setting (29). 

4.3 Option replication with transaction costs in general SVJJ models 

In this subsection we study the problem of option replication presented in Section 3 for general 
SVJJ models (29). We show that in the same hedging policy as in SVJP defined in Section 3, 
jump effects can be ignored in asset as well as in volatility. First, let us recall from Section 3 that 
the hedging error takes the form 

V 1 — h(Sf) = + ^2,n ~ h,n ~ K ^ni 

where I i „, i = 1, 2,3 and T n are defined as in (19). The following conditions on volatility dynamics 
are needed in this section. 

(C 3 ) The coefficient functions ai,i = 1,2 are linearly bounded and locally Lipschizt. Furthermore, 
the common distribution of jump sizes in volatility F v admits the integrability condition 

f z 2 F y (dz) < 00 . 

Jr 


Condition (C 3 ) implies that sup 0<t<1 E y 2 < 00 , which is necessary for approximation procedure. 

Theorem 4.1. Under conditions (Ci) — (Co) — (C 3 ), the limit results in Theorems 3.1 and 3.2 
still hold. 


5 A special case: deterministic volatility with jumps in asset 
price 

In this section we consider a simpler model where asset volatility is a constant. In the absence of 
jumps in asset price, it is well-known that using the classical adjusted volatility 

Sf = a 2 + QoVnf'it), g 0 = cr\/8jn. (30) 
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leads to a non-zero discrepancy between asymptotic portfolio value and the option payoff. In 
particular, [1-1] proved that Vfi converges in probability to 

h(Si) + min(Si, K) — kJ(Si, go), 

where T(a;, g) defined by (22) with cr(y) = a = constant. It was then proven in [25] that asymptotic 
distribution of the normalized corrected hedging error is a mixed Gaussian. In particular, for 
Leland’s strategy and uniform revisions, 

n^ A (Vr - h(Si) - min(5i, K) + kJ{S u g)) 

converges weakly to a centered mixed Gaussian variable for any g > 0. In order to remove the 
corrector in Pergamenslichikov’s result, [5] applied the Lepinette strategy with a 2 defined in (30). 
In fact, they showed that for general European options with payoff function h verifying a power 
decay property, n /3 (P" — h(S\)) converges to a mixed Gaussian variable. Note that these mentioned 
results are obtained in the continuous setting for asset price. 

When jumps are present in asset price, the approximation in these works should be considered 
heavily and seemingly this desired extension is far away from being obvious in their approach. 

One advantage of our method is that it is not difficult to prove that possible jumps in asset 
price can be ignored with very small probability. In fact, we claim that the same result is true for 
jump-diffusion models. 

Theorem 5.1. Assume that the asset dynamics is given by 


N t \ 

+j > 

where N t is a Poisson process with intensity 6 and (£j)j>i is an i.i.d sequence of random variable 
whose common distribution satisfies condition (Ci) and b is a bounded continuous deterministic 
function. Suppose further that the adjusted volatility a is defined by (30) or by the simple form 
(17) - 

Then, for Leland’s strategy 

n /3 ( Vfi - /i(Si) - min (S U K) + kJ(S u go)) 

converges weakly to a centered mixed Gaussian variable. This is still true for any g > 0 in place 
of go and hence the result in [25] is recovered. 

If the Lepinette strategy with a 2 defined by (30) is used then, one gets an asymptotic complete 
replication, i.e. n?{V 1 — h(S\)) converges weakly to a centered mixed Gaussian variable. 

Remark 4. In fact, the following result can be extended to the case when volatility is a bounded and 
smooth deterministic function of t with convex general payoff holding a power decaying condition. 
Hence, the result in [5] can be aslo recovered from Theorem 5.1. 

6 Numerical example 

We present in this section a numerical example using Matlab 2012b. In particular, we assume that 
the asset price follow the jump-diffusion model with stochastic volatility driven by an Orstein- 
Uhlenbeck process dyt = (a — yt)dt + bdWt and volatility function a(y) = yoe v + cr TO j„. This can 
be considered as Hull-White’s model with jumps allowed in asset price. The jump size distribution 
is N( 0,0.2) and jump rate is 9 = 3. Parameters are chosen as £>0 = K = 1, a = —1,6 = 0.2, 
a min = 1, yo = 2, g = yJS/ir. Some paths of the asset price and its volatility are plotted in Figure 
2 . 


d S t = S t - b t dt + odW t (1) 
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Figure 4: Corrected hedging error with different n 


Simulating 500 trajectories by the crude Monte-Carlo, we compute the mean values of the 
corrected hedging error Vj ra — h(S 1 ) — min (Si, K) — kJ(Si, yi, g) with k = 0.001. The convergence 
result is somehow slow, see Figure 4, where mean value of the corrector is 0.2465. Simulation also 
shows that the option price inclusive trading costs converges to the super-hedging price Sq. To get 
a reasonably cheap price we can choose a small value of g but this can leads to a under-hedged 
situation, i.e. the corrector min(Si,A') — g) may be now negative. In case the investor 

accepts risk in his final replication goal the option price can be lowered in the spirit of quantile 
hedging. 


7 Proof of Main Theorems 

As usual, the main results established in three Sections 3, 4 and 5 are just direct consequences of 
some specific types of limit theorem for martingales that we are searching for. For this aim, we 
construct a special approximation procedure following the one in [23]. Our main attempt is to 
prove that jump terms appearing in the approximation can be neglected at the desired rate rfi. 
For convenience, we recall in the first subsection the preliminary setup and refer to [23] for the 
motivation. 


7.1 Preliminary 


Define m 1 = n— 


,(P/A 0 


, 2 / 0 + 1 ) 


and m 2 = n — 


• G*/A 0 ) 


2 / 0 + 1 ) 


where the notation \x] stands 


for the integer part of a number x and l t = In 3 n , I* = In 3 n. Below we focus on the subsequence 
(tj) of trading times and the corresponding sequence (A ■) defined as 


tj = 1 - (1 - j/nY and A= / aldu = A 0 (l - t ) *f >, m 1 <j<m 2 . 


(31) 
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Note that (tj) is an increasing sequence with values in [£*,£*], where i* = 1 — (Z*/A 0 ) 4/3 and 
t* = 1 — (Z*/ A 0 ) 4 ^, whereas (A^) is decreasing in [Z*, Z*]. Therefore, in the sequel we make use the 
notation A tj = tj — tj- 1 whereas AAj = \j~i — A.,-, for mi < j < m 2 to avoid the negative sign in 
discrete sums. 



Figure 5: Two sequences (A j) and (tj) 


Below, Ito integrals will be discretized throughout the following sequences of independent nor¬ 
mal random variables 


^ 1 ,, = 


W t (1) -W t (1) 

3 _ 3 — 1 

~ t j -1 


= 


w m _ w (2) 

3 _ 3 — 1 

\A? ~ t j -1 


(32) 


We set 


p(X,x,y) 


g / ln(a?//v) _ 1\ 
cr(?/) V 2A 4/ 


(33) 


and write for short Pj-i = p(X.;_ 11 S t - ,y t - ). This reduced notation is also frequently applied 

for functions appearing in the approximation procedure. With the sequence of revision times (tj) 
in hand, we consider the centered sequences 


Z 3,j ~ \ Z i,j +Pj- il ~ E (l Z i,j +Pj- il \ Fj-l) ) 
Z 4j = \ z i,j\ - E (\Z hj \ | Fj_j) = \Z hj \ - 


(34) 


The sequences (Z 3 j) and (Z 4 j) will serve in finding the Doob decomposition of considered terms. 
To represent the limit of transaction costs, we introduce the functions 


G(a) = E (\Z + a|) = 2 ip(a) + a (2$(a) - 1), 

A (a) = E (\Z + a\ - E \Z + a\) 2 = 1 + a 2 - G 2 (a), 


(35) 


for a £ K. and Z ~ A/"(0,1). We also write o(n~ r ) for generic sequences of random variables 
(X rt ) n>l satisfying P — lim n _ > , 00 n r X n = 0 while the notation X n = 0(n~ r ) means that n r X n is 
bounded in probability. For approximation analysis, we will make use of the functions 

<t>(\x) =exp| — cj>(X,x) = <j>(X,r}(x)) with rj(x) = \ \n(x/K)\. (36) 
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7.2 Stopping time and technical condition 

We first emphasize that in bounded volatility settings, it is possible to carry out an asymptotic 
analysis based on L 2 estimates as in the previous works [5, 18, 19]. For general stochastic volatility 
frameworks, this approach is no longer valid because k- th moments of the asset prices S may be 
infinite for k > 1 and S is not in general a martingale, see [ 1 , 20]. We come over this difficulty by 
using a truncation technique. In particular, for any L > 0, we consider the stopping time 

t* = t* l = inf jt > 0 : 1 + a t > h} A 1, (37) 

where r\ t = r\{Sf) and a t = max{<r( 7 / t ), |cr'(j/ t )|}. Note that jumps may be not fully controlled for 
stopped process StAr* as in [23]. Therefore, in the presence of jumps we consider its version defined 
by 

St = So exp | b s ds + <r*dW ( s r) - ^ ^ afds + ^ In (1 + Q 1 , (38) 

where of = er (y t ) l{cr(y t )ci,}- Here the dependency on L is dropped for simplicity. Then, it is clear 
that Sf = St on the set {r* = 1}. We easily observe that under condition (C 2 ), 

lim limsupP(r* < 1) = 0. (39) 

L->oo 


For simplicity, in the sequel we use the notation S u = ( S u ,y u ). We carry out an approximation 
procedure for a class of continuously differentiable functions A from R + x R, xR-}K satisfying 
the following technical condition, which is somehow more general than the one proposed in [23]. 
(H) Let A be a R+ x R + xR->l continuously differentiable function having absolutely integrable 
derivative A' with respect to the first argument and for any x > 0, y £ R, 


lim 

n—f oo 

Furthermore, there exist 7 > 


|H(A, x, 2 /)|dA 


C+oo 


|H(A, x, y)|dA =0. 


0 and a positive continuous function U such that 
|H(A, x, y)\ < (A -7 + l)U(x,y)<j>(\,x), 


where cf> defined in (36) and 

sup E U 4 (Sf) < 00 . 

0<t<l 


(40) 

(41) 


Remark 5. In approximation of hedging error, the function U(x,y) takes the form y/x[cia(y) + 
C 20 ’ , (y)] m (up to a multiple constant) for some constants ci,C 2 and m > 0. Therefore, for any 
L > 0, condition (41) is fulfilled as long as sup 0<<<1 ES 1 * 2 < 00 but this is guaranteed by the 
condition of finite second moment of jump sizes (Ci). See Appendix 10. 

For some positive constant L , we introduce the function 

5*0*0 = 5l 0*0 = l a: l 1 {|x|>L- 1 } + ^'iflxlsi- 1 }- ( 42 ) 

Putting rjf = g*(rj t ), one observes that on the set {t* = 1}, 

77 * = L^ 1 and <^(A, S t ) = </>(A, 77 ^) = </)(X, L^ 1 ) := ^l(A), for all t* < u < 1. (43) 
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7.3 Approximation for stochastic integrals 

For the completeness of representation we recall here the asymptotic result established in [23], 
which serves the central role in the proof of the main results. 

Proposition 7.1. Let A(X,x,y) be a function such that A and its first partial derivatives d x A , 
d y A satisfy (H). Then, for i = 1,2, 



Z. j AA J .+ 0 (n-^), 


(44) 


where Aj = A(Xj,S t ~) and A(X,x,y) = A(z,x,y)dz. 

Proof. We follow the argument used in Proposition 7.1 in [23]. Although we are working under 
technical condition (H) which is slightly different from that in [23] but arguments are in fact the 
same. For reader’s convenience let us give the proof in details since the approximation technique 
will be repeatedly used in our analysis. First, making use of the stochastic Fubini theorem one 
gets 

4 = J 4 2 (J A(\,s u )dw^ d t = j ( jT a 2 A(x t , § n )dtj dw«. 

Changing the variables v = X t for the inner integral, we obtain 


/ a 2 A(X t ,S u )dt = A(v,S u )dv = A(X u ,S u )-A(X 0 ,S u ). 

^0 J\ u 

In other words, I n = I l n - / 2 n , where I l n = f* A u dW&, A u = A(X UJ S U ) and I 2n = 
f 1 A(A 0 , S u ) dlfW. Moreover, we have 


h,n = f A u dW.^ + [ * A u dW^ + f A u dW^ := R 1<n + R 2 , n + R 3 , n . 

do jf* j t. 


(45) 


Let use first show that I 2n = o(n /3 ). For any e > 0, one observes that 

P(^l4,nl > 4 < P(^I4,J > e, T* = 1) + P(T* < 1). 

In view of (39), one needs to show that the first probability in the right side converges to 0. Indeed, 
by (H) one has 

/ OO 

e -A/8 dA < CU{x,y)e~ x ° /8 . 

0 

Putting A* = A u and I 2n = J 1 A* dWjfl and making use of the notation S* = ( S*,y ) one 
has 


P (r 


l 2 ,n\ 


> £ ,4 = 1)= P(r 


L 2,n I 


> e) < £ _ 2 n 2/3 E (7*J 2 < C£- 2 n 2 / V A ° /8 sup E U 2 (S*), 


0<t<l 


which converges to zero by condition (H). Hence, I 2 n = o{n P) as n —> oo. Next, let us show that 
R 2 n is the main part of I ln . For this aim, taking into account that l* < X u < A 0 for 0 < u < t*, 
we get R l n = o(n~P). 

Next, let us show the same property for the last term R 3 n in (45). To this end, note again that 
P (n?\R 3 J > e) < P K 1*3,J >s,t* = 1)+P (t* < 1), (46) 
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On the set {r* = 1 } one has the estimate \A U \ < U(S*) f£°(l + z S*)dz = U(S*)f*, where 

f* = /'* (1 + z~ 1 )<j)i,{z)dz. Again on obtains by the Chebychev inequality 

P («^|i? 3>n | > s,t* = 1) = P (n? |^ 3 in | > e) < n 2 ^~ 2 f E(^) 2 d U , 

which is bounded by n 2 ^e~ 2 C{U) ( f*) 2 du with C{U) = sup Q<u<1 Ef7 2 (,$*_) < oo. Taking into 

account that 

/ I pl # / poo \ 2 

Cfnfdu = \- ipJ J (jf (l + z-^M^dz) dX<CX- 4 ^, 

we conclude that lim n _ > . 00 P (n^\R 3 n \ > e, r* = l) = 0 and hence R 3 n = o(n~P) in view of (39). 

It remains to discretize the integral term R 2n using the sequence (Z 1: j). The key steps for 
this aim are the following. First, we represent R 2 n = J),* A u dW^ = Y^=m 1 St 3 i and 

replace the Ito integral in the last sum with Aj_ 1 Zijy / Atj. Next, Lemma 9.1 allows to substitute 
\JAtj = 0 -1 AAj into the last sum to obtain the martingale A im 2 defined by 

k 

M k = g- 1 J2 Aj_iZt jAXj, m 1 < k < m 2 . 

j-m 1 

We need to show that P — lim^oo n^\R 2n — A4 m J = 0 or equivalently, Y^j‘=m 1 B j,n = °( rl_/3 ); 
where Bj n = f* 3 A u jdW^ and A u j = A(X U ,S U ) — A(Xj_ 1 , S t - ). We show this without 
using the Ito’s formula. For this aim, let b > 0 and introduce the set 


tt b 


sup sup 

t* <u< 1 zER 



S u )\+ d x A(z, S u ) 


d y A(z, SJ ) 



Then, for any e > 0, P (n®\ -®j> I > e ) is bounded by 


P(n=)+P(T* <1) + P 


IE 

j=m 1 


B i,n\ 


> £, ty,, T 



Note that lim^^ lim n ^ 0 O P(f2f) = 0 by Lemma 9.4. In view of (39), one needs to prove that 
the last probability converges to zero. To this end, put A u 3 = A u jlnj | <bl5 j and Bj n = 

ft A u ,j dW u ] > where s u,j = \K- | + IS 1 *. - S*_ \ + I y u - - y t -_ \ ■ Then, the last probability 

is equal to P (n^\YS7=m LL | > e) , which is smaller than e~ 2 n 2 ^ Yf7=m E.B 2 by the Chebychev 

\ J 1 J J J J 1 

inequality. Clearly, E B 2 n is bounded by 


2b 2 f’ ((X u - A J _ 1 ) 2 + E (S*. - S) 2 + E {y u - - y t - ) 2 )d u < C ((AA ,) 3 + (A tj ) 2 ) . 

Jtj -1 i- 1 J_1 

Consequently, n 213 Y^'j'Z mi E B 2 n < Cn 2/3 Yl'j= rni (AXj) 3 + (A tj) 2 . Taking into account Lemma 9.1 
we conclude that the latter sum converges to 0 hence, the proof is completed. □ 

Lemma 7.1. Let i(t) = sup{£i : tt < t} and A(X,x,y) is a function satisfying condition (H). 
Then, 
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(i) ■ fo (/‘ (t) K A (K, S u - )du 


dW t (i) = o(n~ 


*), 


* = 1 , 2 , 


(ii). £ (/; (<) A(A„,5 u -)dWi l) ) dW t U) = o(n~ 0 ), i,j G {1,2}. 

Proof. By assumption, |A(A, X, y)| < U(x,y)<p(X,x){l + A -7 ) for some constant 7 and positive 
function U{x,y) verifying (41). Denote by r n the considered double stochastic integral in (i). Put 
A t = /* a 2 A(X u ,S u -)du , we represent r n as 

r„ = [ dA t W t (l \+ f A t dW} x) := r l n + r 2n , 

Jo Jt* 

We will prove that r i n = o(n~ 0 ), i = 1,2. To this end, let L > 0 and consider r* = r£ defined as 
in (37). For i = 1, 2, by r* we mean the ’’corrected” version of r ? - n . i.e. S u , y u should be replaced 
by S* and y* u respectively in A. Now, for any e > 0, 

P (n 0 \ r„| > e) < P (n 0 |r„| > £,r* = l) + P(r* < 1), (47) 

Taking into account X t > l* —> 00 for t £ [0, t*] and using Chebychev’s inequality, one bounds the 
first probability in the right side by 

n 20 £- 2 F,rf n = n 20 £- 2 f EAfdt<Cn 20 s- 2 EU 2 (S*)[ b 2 dt, 

J 0 J 0 

where b t = J* ct ^(1 + A“ 7 )e _A “/ 8 d u. Recall from (17) that 

i[t) 

a 2 u = gy/n(l - = QVn{Xo/XuY, with // = (//- l)/(l + /x). (48) 


Then, splitting the integral as the sum of integrals on the intervals [f J _ 1 , t ;j ] and changing variable 
one gets 


i 20 f b 2 dt<Cn 20 n 2 f d*( 1 + A, u 7 ) 2 e 4 “ d u < Cn 20 3/,2 X(j f d 2 X u 0 (l + A u 7 ) 2 e 


'du, 


•to Jo Jo 

which is smaller than Cn 20 ~ 302 Ag A _ ^(l + A _ 7 ) 2 e _ 4 dA. This implies the convergence to zero 
of the first probability in the right side of (47). In view of (39), one obtains r, n = o(n~ 0 ). Let 
us prove the same property for r 2n . In fact, the singularity at t = 1 requires a more delicate 
treatment. We make use of the stopping time r* again. Put A* = A(X U , S*)!^ i<u(S*)f*}’ 


f* = (1 + A u 7 )^l(A) and r 2 , n = f t l ^ f* a 2 A* u du^J dW^ l \ Then, by the Clrebychev inequality one 
gets P (n 0 \r 2t n\ > e,r* = l) = P ( n 0 |? 2 ,n| > e)- The latter probability is bounded by 


n 20 e~ 2 


sup E U 2 {S* u -) 

0<u<l 



7udu\ dt < Cn 20 e 2 [ ( i(t ) - t) f a*f* 2 dudt := 

J Jf J t (t ) 


On the other hand, for some constant C e independent of n, 


< Cn 20 ~ie~ 2 qXq j' X~?a 2 J: 2 du < C e>( ,n^ J 


A^(l + A u - 7 ) 2 ^ 2 (A)dA, 


which converges to 0 as n —>• 00 . Hence, by taking into account (39) one concludes that P (n 0 |r 2 , n | > e) 
converges to 0. The second equality can be proved by the same way. □ 
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Lemma 7.2. Suppose that A = A(X,x,y) verifies (H). Then, the following asymptotic properties 
hold in probability: 

(i) . £ (/„* 4(A t , SJdW^) d W t ij) = i,j e {1,2}. 

(ii) . f 0 1 A(X t ,S t -)dt = O(n~ 2/3 ). 

(iii) - fo ( ft 1 A(A t ,5 s -)ds) dt = 0 (n- 4/3 ). 

Proof. The procedure used in the proof of Lemma 7.1 can be applied straightforwardly to obtain 
the first equality. Indeed, we can check directly that 




H(A t ,S s -)dW s w ^ cl W t U) = o(n~ 2 ^). 


U) 


Now, consider again the set {r* = 1 } one can prove that (f* A(X t , S*-)dWs^ dW t 

0{n~ 2 P) using again the truncation technique hence, (i) is verified. Next, let us prove (iii) . By 
making use of the change of variable A* = Ao(l — t) 1 /( 4 /3), the double integral is written as 

:= A 0 - 8/3 16/3 2 f X ° A 4/3_1 ( f X ° zW-'Afr^w^dz) dA, v{z) = 1 - . 


By hypothesis, A(X,x,y) is bounded by U(x,y)(l + A 7 )</>(A,ir) for some constant 7 and some 
positive function U verifying (41). Hence, A^Cn is bounded (up to a multiple constant) by the 
double integral 


fV- 1 ^/ A A0 ^" 1 (l + ^ 7 )t/(^( Z /Ao)-)^>^(VAo)-) d "j dA - 

Let oj outside the set {S^- = K}, which has zero probability by Lemma 9.2. It is clear that 
the integrand of the above integral is dominated by a continuous function depending on w, which 
exponentially decreases to 0 at 0 and infinity hence, it is integrable on [0,oo). Therefore, the 
double integral converges to 

J A 4/3_1 17(i§ 1 ) z A P~ x (l + z- J )<j)(z, S' 1 -)d^ dA 

by the dominated convergence theorem. Thus, n 4,s e'„ is bounded in probability. The equality (ii) 
is proved by the same way. □ 


r^o 


7.4 Eliminating jumps 

In this subsection, we establish asymptotic results which will serve in eliminating jump effects in 
our approximation. 

Lemma 7.3. Suppose that 


\B(A,x,y,z)\ < w(z)if(X)U(x,y), for all x > 0, z € R, A > 0, 

where U is a continuous function holding sup 0<t<1 E U 2 (S*_) < 00 for any L > 0 in the definition 
of t* in (37). Suppose furthermore that 


f m 2 (z)v(dz) < 00 and n r f A 4/3 1 (V’ 2 (A) + ip(X))dX —> 0, for any r > 0. (49) 

J r Ji* 
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Then, for any r > 0, 


f [ A(X t ,S t -,y t -,z)J(dtxdz)=o(n r ). 
Jo Jr 


(50) 


Proof. For notation simplicity, one abbreviates B(t, z ) := A( A t , S t -, y t - , 2 ). Let us decompose the 
integral in (50) as 


/ f B(t,z)J (dt x dz) + f f B(t, z) J (df x dz ). 
Jo J r Jt* Jr 


(51) 


Clearly, for any S > 0 and r > 0, P (n r / t » f R B(t, z)J (dt x dz) > <5^ is smaller than P(tVi — 
Nf > 1) = 1 — e~ e A~ l ), which converges to 0. Hence, it suffices to prove the same property 
for the first integral in (51). Indeed, this term can be represented as f* f R B(t,z)J (dt x dz) + 
Jo Jr z)v(dz)dt. We prove that the compensator is almost surely exponentially negligible, i.e. 

n r / B(t, z)v(dz)dt —» 0 a.s. as n —» 00 . (52) 

Jo Jr 

Indeed, by assumption and the change of variable defined in (31), it is estimated by 

Cn r Ag 4 ^ / X 4 ^~ l ip(X)U(S t ^x)-)dX x f w(z)i / (dz), 

Ji* Jr 

which a.s. converges to zero due to (49) and the continuity of U , where t(A) = 1 — (A/A 0 ) 4 ^. Hence, 
it remains to prove that for any r > 0, f 4 f R B(t, z)J (dt x dz) = o(n~ r ) as n —» 00 . To this end, 
note that for any <5 > 0, 


B(t, z)J (dt x dz) 


> 5 < P(t* < 1) + P 


B(t, z)J (dt x dz) 


> 5 , t* = 1 


for any L > 0. Denote by e n the last probability. In view of (39), one needs to show that e ra con¬ 
verges to 0. In fact, by assumption one also has the following estimate | B(t, z)\ < U ( S*_ )ip(\ t )'cv(z) : 
B*(t, z) on the set {r* = 1}. Therefore, applying the isometry for jump integrals yields 


£n = p 

which is bounded by 


B(t, z)1-{\B{t,z)\<B*(t,z)}J x 


>S,t* =1 


n 2r 6 2 E / / B* 2 (t, z)v(dz)dt < n 2r S “ sup E U 2 ‘(S[-) / ip*(Xt)dt x / cc 2 (z)^(dz), 

Jo Jr o<t<i Jo Jr 

and the conclusion follows from (49). □ 


7.5 Limit theorems for approximations 

We first recall the following result in [12], which is extremely useful for studying asymptotic 
distribution of discrete martingales. 

Theorem 7.1. [Theorem 3.2 and Corollary 3.1, p.58 in [12]] Let M n = Yr -1 %-i be a zero-mean, 
square integrable martingale and q be an a.s. finite random variable. Assume that the following 
convergences are satisfied in probability: 

n n 

^e(x 2 1 {| Yi|> 5 } |J- i _ 1 ) —>0 for any 5> 0 and ]T E (X 2 | —> c 2 . 

i—1 i—1 
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Then, the sequence (M n ) converges in law to X whose characteristic function is Eexp(^c 2 t 2 ), 
i.e. X has a Gaussian mixture distribution. 

Below we will establish some special versions of Theorem 7.1. In particular, our aim is to 
study the asymptotic distribution of discrete martingales resulting from approximation (44) in 
Proposition 7.1. 

Let Aj = Ai(X, x,y), i € I := {1, 2,3,4} be functions having property (H) and consider discrete 
martingales (X fc ) mi < fc < m2 and (M k ) mi < k < m2 defined as 

k k 

M k = e~ 1 Yl Y1 A 'j ' Z ' / AA / and = J2 A ij-i z i,j AX j> ( 53 ) 

j—m 1 i£l\{ 4} j=m 1 i£l\{ 3} 

where A i ;] = AfiXj,S t - ) and Zij are defined as in (32) and (34). To describe the limit distributions 
let us introduce 


L = A 2 + 2A 1 A 3 (2$(p) - 1) + A 2 3 A(p) + A 2 2 , L = A\ + A 2 + (1 - 2/tt)A 2 , (54) 

where p is defined in (33). Define now 

2 r+°° _ „ 2 r+°° _ 

C 2 = ftgt^ / A M L(A, S’ 1 -)dA and g 2 = ft gv+i / A M L(A, <S 1 -)dA (55) 
Jo Jo 

with 

M =2 (.T + l)fiTri and ft= (p- l)/(/i + 1). (56) 

Proposition 7.2. Assume that Ai = Ai(X,x,y), i = 1,2,3 and their first partial derivatives d\ Ai, 
d x Ai, d y Ai are functions having property (H). Then, for any fixed g > 0 the sequence (n@M. mo ) 
weakly converges to a mixed Gaussian variable with mean zero and variance g 2 defined as in (55). 
The same property still holds if some (or all) of the functions Ai are replaced by J?° Afiz, x, y)dz. 

Proof. Note that the square integrability property is not guaranteed for the random variables (tv). 
To overcome this issue let us recall the stopping time r* = r£ defined in (37) and put AfiX, x, y) = 
A i (\,x,y)cf>~ 1 (\,x)4>L(X), where </>l(A) defined in (42). Let v* = AfiXp S*_) Z i;j AA^ and 

Ml = E k v*. 

Step 1: We will show throughout Theorem 7.1 that for any L > 0 the martingale M* mn weakly 
converges to a mixed Gaussian variable with mean zero and variance g* 2 (L) defined as 

2 r+°° _ 

g* 2 (L)=ftg / A^L(A, S x - )dA, (57) 

J o 

where L is obtained by replacing all Ai in the formula of L in (54) by the corresponding modified 
functions A i: i = 1,2,3 . To this end, setting a* = E(u* 2 1m „i ,\|J tr ,_i), we first show that 

3 3 s n . ^>o f J 

U o | J 

P (n 2/3 | a *j I > converges to 0. By hypothesis, 


max 

i=l,2,3 


M X u , S*) < U(S*)( 1 + X~ 1 )(j) L (X) < U(S*)(l + A“ 7 ) 


(58) 


for some 7 > 0 and positive function U(S) verifying (41). We observe that 


P j n 2/3 \ ^ a* | > e J = P j n 2/3 \ ^ a*| > e J < e 1 n 2/3 ^ Ea* 


j-m 1 


j-m 1 


J—m 1 
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by Markov’s inequality. Using the Chebychev inequality and then again the Markov inequality one 
gets 



< 9<T 2 (1 + A; 7 ) 4 (AA j ) 4 E U 4 (S* u -)J 2 zf d . 

i= 1 

Taking into account that all of Z i • have bounded moments and using (58) we obtain 

m 2 m 2 

e^n 20 Y Ea * < 9C'£" 1 l r 2 n 2/3 ^ (1 + A,” 7 ) 4 ^) 4 , 

j=m 1 j—m 1 


which converges to 0 by Lemma 9.1. 

Let us verify the limit of the sum of conditional variances E(u* 2 |J r _ 1 ). 

3 j 

A*j_ 1 Z id AXj, one obtains E 
independent. It follows that 


v* .v* .IT - ,- , 
ij 3ji j-r 


) = E ( 


V 2J V 3,j I ^3~ 


■^ = 0 since Z x 


Setting v* . = 
3 and Z 2 • are 


E(u* 2 |J’ i _ 1 ) = E(u* 2 + E(u* 2 l-Fj_i) + E(u 3 * 2 |J- i _ 1 ) + 2 E(u* j u* j |J;_ 1 ). 

Observe that for Z ~ JV(0,1) and some constant a, E(Z |Z + a|) = 2$(a) — 1 and E (Z + a) 2 — 
(E|Z + a|) 2 = A(o), where $ is the standard normal distribution function and A is defined in (35). 
On the other hand, AA ;/ = n _2/3 ( 1 + o(1))/j gi^ X^_ 1 by Lemma 9.1. So, 

n 20 E(u? 2 |.F j _ 1 ) = (1 + o(l))^4r X 0 _ 1 L(A J _ 1 , S*-)AXj. 


Therefore, by Lemma 9.5, the sum n 2/3 V'j 2 E(u* 2 |J r •_-,) converges in probability to ?* 2 (L) 

3—3 J 

defined in (57). Thus, nP weakly converges to A/”(0, <;* 2 (L)) throughout Theorem 7.1. 

Step 2: Let us show that sup e>0 lim^^ limsup n _ >oo P (| n 13 M* m2 — n^M m2 \ > e^j =0. To this 

end, recall that </>(A, S t ) = c/)l{ A) and hence, Ai = A t for i = 1, 2,4 on the set {r* = 1}. Then, the 
conclusion directly follows from 


P (n?\M* m2 - M m2 \ > e) < P (n 0 \M* m2 - M m J > e,r* = l) +P(r* < 1) 


and (39). Moreover, taking into account that <U 2 (L) converges a.s. to c 2 as L —► oo, we conclude 
that n /3 A4 ma converges in law to A/”(0,^ 2 ), which completes the proof. □ 

Let us consider martingales of the following form, resulting from the approximation for Lepinette’s 
strategy, M k = Ej =mi Z i,j + A 2 ,j-i Z 2 ,j + A 4 ,j-i Z 4 j) AA^-. Their limiting variance is 

defined throughout the function 


L(X,x,y) = A 2 (X, x, y) + A 2 (A, x, y) + (1 - 2/tt)A 2 (A, x, y). 


(59) 


The following result is similar to Proposition 7.2. 

Proposition 7.3. Suppose that Ai = Ai(X,x,y), i = 1,2,4 and their first partial derivatives have 
property (H). Then, for any fixed g > 0 the sequence (n^ M. m2 ) weakly converges to a mixed 
Gaussian variable with mean zero and variance g 2 given by (55). The same property still holds if 
some (or all) of the functions Ai are replaced by f?° A)(z,x,y)dz. 
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Proof. The conclusion follows directly from the proof of Proposition 7.2 and the observation that 
E Z 2 . = E(\Z 1 j\ - s/2/n) 2 = 1 — 2 / 7 r, and E (Z i jZ 4 j) = 0, for *=1,2 and m 1 < j < m 2 . □ 

The remaining part of the section is devoted to prove main results following the scheme of 
[23]. Our first step is establish the asymptotic representation at rate vP for each term contributing 
in the hedging error. The approximation procedure also provides the residual parts as discrete 
martingales for which, Propositions 7.2 and 7.3 will be applied to obtain the limit distribution at 
the last step. 


7.6 Approximation for I ln 

The following approximation is obtained in [23]. 

Proposition 7.4. Let H = (z -1 / 2 /2 — z~ 3 / 2 ln(x/K))tp(z,x)dz and define 


U 


.,k = Q 1 H CT (2/t 3 _i ) s t-_ 1 Hj-i z i,j AA r m i<J< m 2- 


3=m 1 


Then, under (C x ) and (C 2 ), P — lim n _ yoo |/ 1 n — 2min(5' 1 , K) — U lm2 | =0. 

Proof. By (19), one represents I ln as 

h,n = [' d 2 S 2 .C xx {t,S t -)dt- (\ 2 {y t )S 2 .C xx {t,S t -)dt. 

Jq Jo 

The last term is nP negligible by (ii) of Lemma 7.2. To study the first integral let us introduce the 
function v4(A,a:) = x 2 C xx {t,x) and split it as 

f 1 a 2 S 2 -C xx (t, S t -)dt = I' d 2 Sl-C xx {t, S x -)dt + f 1 a 2 (A(X t , S t ) - A(X t , SJ) dt. 


The first integral f a 2 S 2 _C xx (t, S 1 -)dt almost surely converges to 2min(S' 1 -, K) faster than n r 
for any r > 0, see [23]. Let us study the last term which describe jumps of A. Using the Ito Lemma 
for A( A t) S t ) — A(Xt, S 1 ), we rewrite it as 


,„+ / 1 I' d 2 d x A{\ t ,S u -)a{y u )S u -dW^dt, 
Jq Jt 


(60) 


where 


1 pi 


a 2 A 1 (\ t ,S u -,y u -)dudt and 


0 


A(Xt,S u -, z)J(du x dz)df 


with 


Ai{\,x,y) = d t A(X,x) + d xx A(X,x)a 2 (y)x 2 , A(X,x,z ) = A(A,a;(l + z)) - A{X,x). 

Then, the approximation procedure of Proposition 7.1 is used to get a discrete martingale approx¬ 
imation U 1 m2 for the Ito’s integral of (60). 

Now, let us show that e i n = = 1,2. In fact, e ln = o(n _/3 ) by (iii) of Lemma 7.2. The 

jump term e 2 n can be represented as e 2 n = J 1 J R ' a 2 A(Xt, S u - , z)dt) J(du x dz) by Fubini’s 
theorem [2]. Changing variable v = d 2 dt as in (31), one gets 

f a 2 A(X t ,S u -,z)dt = [ A(v,S u -,z)dv := D(X u ,S u -,z) 

Jo Jx u 
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and hence, e 2 n = J 1 J R D(X u , S u -, z)J(du x dz). On the other hand, 


r^o 


pS — (I+ 2 ) /*A 0 


^(Au, ,z)= A(X U , S u -, 0 ) = 


/s _ 


<9;,A(«, x)dt>dx. 


Direct computation shows that d x A(v, x) = 2xC xx (v, x) + x 2 C xxx (v, x) and 


1 


1 


C*x(^ x ) = —■/=¥>(«> x )> C xxx {v,x) = ~ — <p(v,x) -y/v + 


In (x./K) 


Denoting ip(v,x.) = -^(f> 0 (v)e ^ with <f>o{v) = y v / 8 and using the fact that v k e v / 2 

is uniformly bounded for all k, one has \d x A{v, x)| < C-k=( l+u _ 1 )</>o(u), for some positive constant 

V x 

C. This estimate implies that 


r^u~ ( 1 +*) \ 

Is _ 7^ <lx 


r*A 0 


x / (1 + u Ofo^dv <Cm(z)</>o{X u )S u -, (61) 


I D{X u ,S u -,z)\<C 
where 

1 ~ f°° 

m(z) = zl {z>0} + | r m /2 1 {-i<g<o} + !> ^°^ = J X ( 1 + v )M v ) dv - (62) 

Clearly, </> 0 and to satisfy condition (49) of Lemma 7.3 hence, e 2 n = o(n~ r ) for any r > 0. □ 

7.7 Approximation for I 2n 

Proposition 7.5. Under (Ci) and (C 2 ), ri^I 2n converges to 0 in probability as n —► 00 . 

Proof. We represent I 2 n as 


f 1 a(y t )S t -A(t)dW^ + [ [ zS t -A(t~)J (dt x dz ) := b 1>n + b 2>n , 
J 0 Jn JR 


(63) 


where Aft) = C x (L(t),S^ t \) — C x (t, S t ). We first claim that the Ito’s integral of (63) can be omitted 
by Lemma 7.1. To see this, it suffices to apply the Ito’s formula, one represents the difference A t 
as 

/' (d xt (u,S u -) + a 2 (y u )S 2 u -C xxx (u,S u -))du+ f C xx (u, S u - )a(y u )S u -dW^ 

A(i) J L {t) 

+ [ [ (d x (u,S u -0 + z))-C x (u,S u -))J(dzxdu) 

J i(t) JVL 


In view of (16), 


C xt {u,x) = --d 2 u (2xC xx (u,x) + x 2 C xxx (u,x) S ) := d 2 u A(u,x). 


(64) 


Therefore, b, n equals the following sum 

[ [ vlv(yt)S t -A(u,S u -)dudW t {1) + f f a{y t )S t -a 2 (y u )Sl-d xxx (u,S u -)dudWj ) 

Jo J i(i) Jo J i(t) 

+ f f [ *(yt)S t -(d x (u,S u -(l + z))^d x (u,S u -))J(dzxdu)dW t W . (65) 

•to Jdd 
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The first two integrals converge to 0 more rapidly than n ^ by Lemma 7.1. Let us study the jump 
term in (65), which will be denoted by a„. Clearly, by the Fubini’s theorem a n equals 


E 

l<i<r 


'&(u,S u -,z) 


a(yt)S t -dW ^ ) J(dz x du) 


( 66 ) 


where ^(u, x , z) := C x (u , a;(l + z)) — C x (u , a;). We prove that a„ = o{n~ r ) for any r > 0 following 
the demonstration of Lemma 7.3 with some modification. In particular, we decompose the sum in 
(66) into two parts: a x n , the first concerns the index i with m 2 < * < n and the second one a 2 ra , 
which is the sum over the rest of index i, 1 < i < m 2 . Clearly, P(n r ’|a 1 „| < S) < P (Ni — N t . < 
1) = 1 — e -0(i-* \ which converges to 0. 

To study the second one a 2 n , we run again the argument used to to obtain the estimate (61). 
In particular, z)\ is bounded by 4>i(A u ) f x < 0 1 (A u )y / xzuo(z), where 


</>i(A) = y/K/{ 27r)A t 1//2 e A ^ 8 and vj 0 (z) = 1 + 


Denote by a 2 the compensator of a 2 n . Then, it is clear that 


\/T" 


L { — 1 < 2 < 0 } 


(67) 


*U< E / Mn,S u -,z) 

l<i<m2 b-l 

< e r 

1 <i<m 2 Jt i -1 


a(y t )S t -dW t W 


v(dz)du 


a(y t )S t -dW t W 


du x -coQ(z)v(dz) 


It is important to highlight that X u := 


( 68 ) 


may not be squared integrable. 


«- \C r(y t )s t -dW t {1) 

To overcome this issue, consider the stopping time r* dehned in (37) for some L > 0. On the set 
{t* = 1}, one has for u £ 


E X 2 U * = E .fsX / ’ v(y* t )S* dW t (1) ) < ES*_ / ‘ S**dt < CL 


-2-1 


and hence, EX* < i/EX^* < CLn 1 / 2 by Cauchy-Shwart’s inequality. Therefore, 


P(n r K,„l >8,T* = l)<n r S~ i ]T 


<^i(A u )EX*du x / m 0 {z)v(dz) 


l<i<rri2 1 


< n r 5 1 CLn 1 ^ 2 / <(>i(A u )d u x / ^(-^(dz) 

l<i<m 2 "El , ' R 

< n r <W 1 C'Ln“ 1 / 2 f cj)i(A u )du x f ujQ(z)v(dz) 

J 0 J R 


(69) 


Now, taking into account -a7o(^)^(dz) < 00 and J* ^>i(A u )du rapidly converges to 0, one con¬ 
cludes that the right side of (68) converges to 0 for any r > 0 and so is P(n r |a 2 J > S,T* = 1). 
Noting that 

Pfa’KJ > 6 ) < PK|a^| > <5,r* = 1) + P(r* < 1) 

and using (39) one obtains n r a 2 —>• 0 in probability for any r > 0. 

Now, putting u 2n = a 2n ~ a 2 n > we nee d t° show that P(n r jS 2n | > S) —>■ 0. To this end, 
consider again the stopping time r* defined in (37) for some L > 0. On the set {r* = 1}, 
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one has \t&(uS u -x, z)\ < S^_<fi(X u )wo(z), where S'*. is the stopped version of S u ~. Clearly, 
sup 1<i<n sup t < u<t . E S*_ | Wf 1 ^ — WW| 2 < Cn _1 for some positive constant C. It then follows 
by the Chebychev inequality that P(n r |a 2 n \ > 5, t* = 1) < n 2r S~ 2 E5* 2 n , where 5 2 2 n is obtained 
by substituting S u - by S*_ in the function d>(u, S u -, z). Now, the well-known isometry for jump 
integrals applying to 5 2 n = a 2 n — implies that E a* 2 n is bounded by 


E E 

l<i<m2 


\*{u,S*,z)\ 2 \WW - wW\ 2 v(dz)du, 


which is smaller than 

rt 


f ^i(A„)EX* 2 d u x f WQ(z)iy(dz) < Cn 1 f (ff(X u )du x f n7g(z)i/(d^:). 

J o Jq «/R 

Again, f R zcro(z)is(dz) < oo by condition (Ci). Therefore, for some constant depending on L, 

P(n r |5 2 r! | > S, t* = 1) < C(L)n 2r S -2 ^ 1 f (f> 2 (X u )du x f zu 2 (z)iy(dz), 

J o «/K 

which converges to 0 for any r > 0. Letting now L —> oo and using (39) we obtain that \a 2 n \ = 
o(n~ r ) for any r > 0. By the same way, one can show that n r b 2 n —> 0 in probability for any r > 0 
and the proof is completed. □ 


7.8 Approximation for J 3n 

Proposition 7.6. Suppose that (Ci) and (C 2 ) hold. Then, for any r > 0, n r \I 3n \ —> 0 in 
probability as n —> oo. 

Proof. By (19), one has B (t,S t -,z) = f$ C xx (t,u)dudv. Recall that C xx (t,u) = 

u~ 1 Xf 1 ^ 2 tp(X t , u) < u~ 3 ^ 2 cj)i(Xt), where 4>i(X) = 2 r K)Xf 1 ^ 2 . Direct calculus leads to 

1 /2 

|B(f,S t -,z)| < CS t L cfi(Xt)w(z) where w(z) defined in (62). Therefore all assumptions in Lemma 
7.3 are fulfilled and the conclusion follows. □ 


7.9 Approximation for r n 

Let us study the trading volume L n . It is easy to check that for v > 0, 1 — $(u) < Cv^ipiv) and 
Jo p(X u , S u )du + f t <p(X u , S u )du almost surely converges to 0 more rapidly than any power of n. 
Therefore, one can truncate the sum and keep only the part corresponding to index ni\ < j < m 2 . 
Next, one can ignore jumps terms that may appear in approximations via Ito’s formulas in the 
interval [t *, 1]. For convenience, let us recall here the approximation result for T n obtained in [23] . 

Proposition 7.7. Under conditions (Ci) — (C 2 ), the total trading volume T„ admits the following 
asymptotic form 

r n = r(Si, 2 / 1 , e) + (^2,m 2 +^3,m 2 ) + o(« _/3 )- 

7.10 Proof of Theorem 3.1 

By Propositions 7.4-7.7, the hedging error is represented as V{ l —h(S i) = min (Si -, K)—kT (S 1 -, y±, p)+ 
A4 m2 , where the martingale part of the hedging error is given by M k = \U lk — k(U. 2 k +U :i k ) and 
hence, the sequence (n 13 A4 m2 ) converges in law to a mixed Gaussian variable by Proposition 7.2 
and Theorem 3.1 is proved. □ 
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7.11 Proof of Theorem 3.2 


Suppose now that the Lepinette strategy 7 " is applied for the replication problem. In the same 
principle one can represent the corresponding hedging error as V™ — = \l ± n +I 2 n ~l 3 „ —fiTn, 

where 

rU-i ^ 

L 2,n = h,n + AS t, / C xt( U i S u) du 

i> 1 J o 

and f„ = 1 \lt I i s the trading volume. Recall that I 2 n est negligible by Proposition 

7.5. Let us investigate the above sum. By (64), it can be represented as 

[ Mu,S u )dv j a(y t -)S t -dW t {1) f A(u,S u )dvj zdS t ~J(dz x dt) 

€>1 ^1 Jti —1 i>l J ti~ 1 

using the usual change of variable, where A defined in (64). Now, the approximation technique of 
Proposition 7.1 can be applied to replace the first sum by martingale U 2 m2 defined by 

k 

u 2 ,k = 8 _1 Y j -1 Z i,j AX j > rni<k<m 2 

j=m 1 

and Y (A, a;) = z -3 / 2 In (x/K)ip(z, x)dz. On the other hand, one obtains the same estimate (61) 

for the integrand, which implies that the second sum can be omitted at order n r for any r > 0 by 
Lemma 7.3. 

Now, approximation representation for the trading volume P ra following the procedure in the 
approximation of r„. The following is established in [23]. 

Proposition 7.8. Under conditions (Ci) — (C 2 ), 

P - lim n p \T n - 77 min (SI, K) - (U 2 +U 3 )\ = 0. 

n—> 00 

Hence, = IA { + U l m — K*(U 2t m 2 +^ 3 ,m 2 ) is the martingale part of the hedging error for 

Lepinette’s strategy, which can be represented in the form 

k 

■Mk = 8 1 ^2 {A l j _ 1 Z 1 j + A 4j Z 4j _ 1 + A 2 j_ 1 Z 2 j)AAj 

j=rn 1 

for explicit functions A i holding the assumption of Proposition 7.3. Then, the convergence in 
law to a mixed Gaussian variable of the sequence (n 13 Ai m ) is guaranteed by Proposition 7.3 and 
hence, Theorem 3.2 is proved. □ 

7.12 Proof of Theorem 4.1 

Note first that the approximation representation for the replication error is the same as in SVJ case, 
in particular, approximations of 1^ n , i = 1 , 2 ,3 are the same since martingale sums are resulted 
from Ito’s formula in one dimension case. The only difference is that in finding the limit of the 
total transaction costs one has replaced S t - and y t - by terminal values S 4 - and y 4 ~- Now, the 
two-dimension version of Ito’s formula applied for the difference provides sums concerning to the 
dynamics of y*. By the elementary property of Poisson process one can ignore the jump part of 
y t in the time interval [t*, 1]. Hence, the martingale approximation for this difference is the same 
as in SVJ case. However one needs to check the integrability of ai(t,y t ),i = 1,2. For this aim, 
condition sup 0<t<1 Ey 2 < 00 is needed but this is fulfilled under condition (C 3 ) together with the 
linear growth and Lipshitz properties of these coefficients, see Appendix 11. □ 
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7.13 Proof of Theorem 5.1 


The proof can be proceeded by the same method as that for Theorem 3.1 and Theorem 3.2 but 
with more simple arguments. In fact, the difference between the use of a 2 = a 2 + goy/ n/'(f) and 
that of the simple form a 2 = go yj nf (t) only comes from the approximation due to the substitution 
A = J ) 1 a 2 ds. In particular, 

X t = a 2 (1 - t) + gy/n y/f'(t) = cr 2 (l -t) + A t , 

where A* is defined by the same formula with the simple form. Then, Lemma 9.1 is still true for 
the sequence (A_,) constructed from the classical form (30). Of course the index mi, m 2 are now 
replaced by mi, m 2 defined by 

m 1 = n — [ra9 _1 (r)], m 2 = n — [ra? _ 1 (Z*)], 

where ’d(z) = a 2 z^ + Ao 2 *- 1+M ^ 2 , which is an increasing function for fj > 1 , and the notation [x] 
stands for the integer part of a number x and l* = In -3 n, l* = In 3 n. Similarly, we consider a 

subsequence (tj ) of trading times and the corresponding sequence ^A ? j defined as 

tj = 1 - (1 - j/nY and X j = cr 2 (l - tj) + A 0 (l - > fh 1 < j < m 2 . (70) 

The rest of approximation procedure is the same as that of Theorem 3.1 and Theorem 3.2. □ 

8 Concluding remark 

Diffusion-based stochastic volatility models well account for volatility clustering, dependence in 
increments and long term smiles and skews but can not generate jumps nor realistic short-term 
implied volatility patterns. These shortcomings can be fixed by adding jumps into the model. There 
are two possible ways to emerge jumps into stochastic volatility models: adding an independent 
jump component to the return or in the volatility process itself. We showed that jumps in such 
frameworks do not affect asymptotic property of the replication error in approximate hedging 
with transaction costs. The results established in the present note is general enough for practical 
purposes. 

It should be mentioned that in [29, 30], the authors studied the asymptotic property of hedging 
error resulting from discrete delta hedging in exponential Levy models without transaction costs. 
More precisely, they showed that the normalized hedging error converges stably in finite-dimension 
laws to an explicit variable. This result was applied to a problem of hedging a discontinuous payoff 
option in Merton’s jump-diffusion model. However, they left the jump residual as an unhedgeable 
term even in sense of approximate hedging. 

It would be interesting to investigate the problem of approximate hedging options written on 
multiple assets where jumps are allowed in asset price processes. In the absence of jumps, such 
a problem has investigated in [24] and hence, it is reasonable to believe that jumps influence can 
be also removed in the limiting property of the hedging error in multiple frameworks. Another 
interesting problem is to study asymptotic properties of jump risk in small transaction costs models. 
In fact, when k = K 0 n~ a for 0 < a < 1/2 the complete replication can be obtained for both Leland 
and Lepinette strategies for deterministic or local volatility models. Such extensions are in progress 
research. 
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9 Auxiliary Lemmas 

Lemma 9.1. There exist two positive constants C\. C? such that 

v 0 {h) < inf | AAj-1 < sup \A\j\ < C^toT 2 ^vq(1*), (71) 

m i—i — m 2 777.^ <j<rn 2 

where vq{x) = Moreover, let jl = \(p + l)/x' r + T , then 

AXj = /ln _ 2 / 3 i/o(Aj_i)(l + o(l)) and AXj (At ,-) -1 / 2 = g(l + o(l)) (72) 

Proof. It follows directly from relation (31). □ 

Lemma 9.2. For any K > 0 and 0 < t < 1, P (St = K ) = 0. 

Proof. We prove that for 0 < t < 1 and any real number a, P('i/*t = a) = 0, where ift = 
J q b t d t + f* a (y s ) dlT'i 1 ' 1 - \ f* a 2 (y s ) d s + 1 R (1 + £j) ■ Indeed, one can represent = 

pB t + y/l — p 2 Z t , where B t is the Brownian driving y t and Z is another Brownian independent 
of B. Now, conditionally on the Brownian B and jump terms J^^iln(l + £j), i/t is a Gaussian 
variable. □ 

Lemma 9.3. For any s > 0 and K > 0, limsup : P (inf 1)<t<1 1 In (S t /K)\ < s) = 0. 

Proof. Let p be some positive number. Clearly, the above probability is bounded by 


P( inf \ln(S t /K)\<e,N 1 -N v = 0,\ln(S 1 /K)\>p) 

V<U< 1 

+P(| ln(jS' 1 /^sT)| < p) + Pim -N v > 1 ). 


(73) 


Let us show that the probability in (73) is equal to zero for v sufficiently close to 1. On the set 
{Ni — N v = 0}, we have In (S t /K) = ln(S 1 /K ) — if t , where ip t = —^ f * a 2 (s)ds + /J er(s)dW,,(s). 
We can check directly that tp*(v) = sup < <1 \if r {u)\ —> 0 a.s. as v —> 1. So, if |ln(S’ 1 //\)| > p 
then for v sufficiently close to 1 and v < u < 1 , 

|hW^)| = MSi/K) - A\ >\\HSi/K)\ - Wl > \ \HSi/K)\ > p/2. 


Therefore, for p > 2e, one obtains inf t , <ll<1 | ln(5 t /iir)| > p/2 > £ and so, the first probability in 
(73) is equal to 0. On the other hand, P(W — N v > 1) = 1 — e _e ( 1 ~' l d —» 0 as v —» 1. Letting now 
p —> 0 we get P(| ln(S' 1 //\)| < p) —>• P(5 1 = K ) = 0 in view of Lemma 9.2 and hence Lemma 9.3 
is proved. □ 


Lemma 9.4. Suppose that A = A(A, x, y) and its partial derivatives d\A, d x A , d y A verify condition 
(H). Set 


r n = sup (\d x A(z,r,d)\ + \d x A(z, r,d)\ + \dA(z,r,d )\), 

(z,r,d)£ [l* ,/*] xt3 


where B = [S min -, 5 max -] x [y mi 


min' ’ 2/max 


with 


S min~ = .. i nf ,. 5 max^ = SU P S u~ > S/min" = .. K. Vu~ - 


t* <u<£* 


t* 


t* <u<£* 


= SU P Vu~- 

t* <u<t* 


Then, lim^^ lim„^ 00 P(r n > b) = 0. 

Proof. See Lemma A.4 in [23] with remark that the left continuity of S t - and y t - gives the same 
argument. □ 
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Lemma 9.5. Suppose that A = A(A,x,y) and its first partial derivatives have property (H). Set 
A(A, x, y) = f x A(z, x, y)dz and A(A, x,y) = A (A, x, y). Then, for any 7 > 0, 


P - lim 

n—> 00 


' '3 

j=m 1 


A7 1 A(A j _ 1 ,S t - )AA j- A 7 A(A, 5' 1 -)dA 


= 0 , 


where S t = ( S t ,y t ). The same property still holds if A(A,x,y) = A(A,x,y) or the product of these 
above kinds. 

Proof. See Lemme A.5 in [23]. □ 


10 Some moment estimates 

Lemma 10 . 1 . Let yt is some ltd’s process and St be the asset process given by 

St = So exp | J b t dt + J a {y s )dW ( s 1) J a 2 (y s ) dsj (1 + £,), 

where N t is a homogeneous Poisson process with intensity 9 independent of a sequence of 

i.i.d. variables. We assume that the jumps ingredient (£j)j> 1 , N t are independent of the Brownian 
motion W t and of that of y t . If b and 0 are two bounded functions then, for any m > 0, 

E57 < C(m) exp{0£(E(l + £i) m - 1)}, for all t G [0,1], 

where C(m ) is some constant depending on m. 

Proof. Let us represent S t = b t £t(cr)X t with X t = nj=i (1 + £j) and 

b t = S 0 eti b ° ds , f t (cr) = exp |y a (y s ) dVF s (1) - ^ ^ ct 2 ds } ’ 

By hypothesis, the stochastic exponential £ t (er) is a martingale with expectation 1, indepedent of 
X t . Therefore, ES™ < CE£f l (a)EXfi 1 since sup 0<t<1 6 ™ < C. Because a is bounded one has 

E £f n (cr) = E£ t (mcj)e( m2 - m)/2 ^ 2 (^) ds < CE£ t (ma) = C. 

On the other hand, using the usual conditioning technique gives 

N t 

EA7 = E H (1 + = exp{0f(E(l + &)"* - 1)}, 

3 =1 

which implies the desired conclusion. □ 

Lemma 10 . 2 . Under the assumptions of Lemma 10.1, for 0 < u < v < 1, 

E(S U - S v ) 2 < C\u - v\, 


for some constant C. 
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Proof. For 0 < u < v < 1, put b v/u = b ‘, X u/v = II^jv„+i ( x + 0) and 


£„/Jcr) =exp {/ <T (y°) dW s 1) -7} cr 2 ( 2 / s )ds|. 

Then, £ v / u (cr) and are independent and 

sup (E£* u (<j) + Ee 2d “ b ”) < oo (74) 

0<w<i><1 

since b ans a are bounded. Denote 5 = 0(v — u). It is easy to check that 

E(X u/v - l ) 2 = e 5(ESl+E{ ^ - 2e 5E£l + 1. (75) 

Let us first show that Ei{X u j v — l) 2 < CS , for some constant C. Obviously, for any finite interval 
[a, 6], \e x — 1| < Cx by Taylor’s approximation. From condition (C 3 ), E£ 2 < oo. Now, if E£i = 0 
then E(X u /„ — l) 2 = e 5E£ i — 1 < CS. Similarly, in case E£- ( + E£ 2 = 0 one has E£ 3 ^ 0 and hence 
E(X u/v -l) 2 =e^-l < CS. Lastly, if both of E£i and E£i + E£ 2 are non zero one can estimate 
E(X U / V — l) 2 by |e 5 l E£l+E£ i — 1| + 2|e <5E£l — 1| < CS. Using the same argument one can easily 
prove that 

E(4/» - l ) 2 < C6 and V(b v/u - l) 2 < CS 2 . (76) 

Clearly, E (S u - S v ) 2 = E5 2 E - l ) 2 and 

(f; - l) <2 (& 2 /„(4/» - l ) 2 + (b v/ u - l ) 2 +bliA lv (<T)( X u /v ~ l) 2 ) ■ (77) 

By Lemma 10.1, sup 0<u<1 ES 12 < oo. Now, taking expectation in (77) and using (74), (75) and 
(76) one obtains the conclusion. □ 

11 Stochastic differential equations with jumps 

In this section we recall the basic result in the theory of stochastic differential equations with jumps 
(SDEJ) of the form 

d y t = (£, yt)dt + a 2 (t, y t )dW t + d£t, (78) 

on the time interval [0, T] with initial value yo , where Q = £j is a compound Poisson process 
independent of the Brownian motion W and Ej/q < oo. 

Theorem 11.1. Suppose that on, i = 1,2 are locally Lipshitz and linearly bounded functions and 
Ez/g < oo. Assume further that jump sizes of Qt have finite second moment. Then, there exists a 
unique solution y t to (78) with initial value y$ and 

E ( sup yt) 2 < C(T)(1 + Ez/g) < oo. (79) 

0 <t<T 

Furthermore, for any 0 < s < t < T, there exists a positive constant C such that 

E|z/ t -z/ s | 2 <C|t-s|. (80) 

Proof. The existence and uniqueness of the solution follows by adapting the classical method used 
for SDEs, see for instance Theorem 2.2 in [10]. To prove (80), we note that 

/ * \ 2 

Ea 2 (M, 2 /«)du + 3E ^ I 

\j=N s +1 ) 


E|y t - y s \ 2 < 3E / ai(u,y u )du) +3 
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By the linear boundedness of aq, a 2 and (79) one gets 


E 


oil {u, y u ) d u^j < C\t - s\ (1 + Ey 2 )dw < C\t - s[ 


Similarly, Ea| (u, y u ) d u < C f*( 1 + Ey 2 )dit < C\t — s|. To compute E +1 €j) we apply 

the conditioning technique to get E ^ ^j— N +1 — s\Var(;i + (E£i) 2 (A|t — s| + A 2 |t — s| 2 ) 

and the conclusion follows. □ 
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